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Motivated by the experimentally observed high mobility of gold atoms on graphene and their 
tendency to form nanometer-sized clusters, we present a density functional theory study of the 
ground state structures of small gold clusters on graphene, their mobility and clustering. Our 
detailed analysis of the electronic structures identifies the opportunity to form strong gold-gold 
bonds and the graphene mediated interaction of the pre-adsorbed fragments as the driving forces 
behind gold's tendency to aggregate on graphene. While clusters containing up to three gold atoms 
have one unambiguous ground state structure, both gas phase isomers of a cluster with four gold 
atoms can be found on graphene. In the gas phase the diamond shaped Auf cluster is the ground 
state structure, whereas the Y shaped Auf becomes the actual ground state when adsorbed on 
graphene. As we show, both clusters can be produced on graphene by two distinct clustering 
processes. We also studied in detail the stepwise formation of a gold dimer out of two pre-adsorbed 
adatoms, as well as the formation of AU3. All reactions are exothermic and no further activation 
barriers, apart from the diffusion barriers, were found. The diffusion barriers of all studied clusters 
range from 4 to 36 meV, only and are substantially exceeded by the adsorption energies of -0.1 to 
-0.59 eV. This explains the high mobility of Am -4 on graphene along the C-C bonds. 

PACS numbers: 71.15.Nc, 78.67.Bf,68.43.Fg 



I. INTRODUCTION 

Historically, carbon allotropes, e.g. graphite, carbon 
nanotubes, and fullerene have been extensively used for 
studying the absorption processes of finite-sized particles 
on them. Small coinage metal clusters on single layers 
of carbon, i.e. graphene, were first used to model the 
adsorption on graphite surfaces or single-walled carbon 
nanotubesP^ After the experimental evidence for the ex- 
istence of graphene was found by Novoselov et al. 4 5 , now 
even awarded with the Nobel Prize in Physics, the study 
of cluster adsorption on graphene became important in 
its own right. 

Graphene stands as an extraordinary material that of- 
fers enormous possibilities for applications in electronics, 
sensors, biodevices, catalysis, and energy storage P*"^ The 
unique electronic structure of graphene with a linear dis- 
persion of the electronic structure at Dirac points not 
only plays an important role in charge transport but in 
verifying theoretical predictions in quantum electrody- 
namics by table-top experiments, too. As graphene is a 
zero band gap material, one of the routes to make use of 
graphene in electronics industry is to create a small band 
gap by means of functionalization with external chemical 
agentsPHIU 

By means of density functional theory, the effect of Au 
nanoparticles on the electronic structure of graphene has 
been studied.^ The authors showed that AU38 nanopar- 
ticles covered with methylthiolate molecules introduce 
new Dirac-type points due to charge transfer but keeping 
the graphene layer metallic whereas the bare nanopar- 
ticle opens up a small band gap at the Dirac point of 
graphene. A recent workP^ has emphasized the effects of 
the distribution of Au atoms on the charge carrier mobil- 
ity of graphene. They have concluded that the formation 



of Au clusters increases the mobility whereas a homoge- 
neous distribution reduces the mobility. 

The adsorption and diffusion of gold adatoms on highly 
oriented pyrolytic graphite (HOPG) have been studied 
for decades, first, by means of the desorption flux from 
the surface 15 , later on directly with scanning tunneling 
microscopy 16 , and transmission electron microscopy™ 
Also ab-initio calculations on the adsorption of small 
gold clusters and the diffusion of adatoms and dimers 
on HOPG have been reportedP^o] 

Very recently the growth of gold nanoparticles on few 
layers of graphene by physical vapor deposition has been 
monitored. 21 The growth process stopped at an average 
particle size of 6.46±0.68 nm on a single graphene sheet. 
As atomic-scale manipulation is nowadays possible in ex- 
periments, the study of small clusters and their effects on 
the electronic properties is highly relevant for the future 
realizations of nanoscale devices. 

A systematic theoretical study of the size variation 
of small Au clusters adsorbed on graphene is missing. 
Here we address this issue by means of density func- 
tional theory (DFT). We study the adsorption of Aui_4 
on graphene, focusing on an analysis of the electronic 
structures. Diffusion barriers are calculated from the to- 
tal energies and confirm the high mobility of Aui_4. We 
study in detail the step-wise formation of a gold dimer 
out of two pre-adsorbed adatoms, as well as the forma- 
tion of Au3 5 4. We show how both AU4 isomers from the 
gas phase can be produced on graphene by two distinct 
clustering processes. 

This paper is organized as follows. In the next sec- 
tion |n| we summarize the computational details of this 
work. In Sec. Ill the ground state structures of 



Aui_4 /graphene, their adsorption energies, charge trans- 
fers, and electronic structures are discussed. Starting 
from these structures the mobility, i.e. diffusion barri- 



2 



ers, and clustering processes have been studied and are 
presented in Sec. [TV) Finally, Sec. |V| summarizes and 
concludes the paper. 



II. COMPUTATIONAL DETAILS 



The scalar-relativistic ab-initio DFT calculations 
were perfor med using the projector augmented wave 
(PAW^PEH method as implemented in vaspP^ The 
exchange-correlation interaction was treated in the gen- 
eralized gradient approximation (GGA) in the parame- 
terization of Perdew, Burke, and Ernzerhof (PBE)!^ A 
cut-off energy of 600 eV was used and a Gaussian smear- 
ing with a width of a = 0.05 eV for the occupation of 
the electronic levels. Spin-polarization was taken into 
account for all calculations. 

The graphene sheet was modeled by a 5 x 5 supercell, 
i.e. 50 carbon atoms, using the calculated C-C bond 
length of 1.42 A. The repeated sheets were separated 
from each other by 20 A of vacuum. A Monkhorst-Pack 
T-centered 5x5x1 k-point mesh (13 k-points in the ir- 
reducible wedge of the Brillouin-Zone) was used for the 
structural relaxations. The relaxation cycle was stopped 
when the Hellman-Feynman forces had become smaller 
than 5 • 10 -3 eV/A. We also used a finer k-point mesh, 
i.e. 16 x 16 x 1 k-point mesh (130 k-points) and found 
the changes in the geometry of the systems to be negli- 
gible. Therefore we decided to use the relaxed structures 
obtained with the 5x5x1 k-point mesh for total energy 
(Eq) and density of states (DOS) calculations. However 
we used a finer k-point mesh, i.e. 20 x 20 x 1 with 202 k- 
points, to accurately calculate E and the DOS of these 
structures. Note that cluster adsorption energies, £? a d s = 
£?o(Aui_4 /graphene) — £?o(Aui_4) — (graphene), are 
negative when the adsorption is exothermic. The pro- 
jected DOS were calculated within Wigner-Seitz spheres 
of radii 0.863 A (carbon) and 1.503 A (gold), using for 
illustrative purposes a higher Gaussian smearing of 0.1 
eV. 

In order to obtain the diffusion barriers the total en- 
ergies of the gold clusters on different binding sites were 
calculated. In the calculations the x-y coordinates of the 
gold atom binding to carbon, as well as the carbon atoms 
in the rim, were fixed, while the rest of the structure 
could fully relax. Since the calculated diffusion barriers 
are less than 36 meV, we checked their dependence on the 
used GGA exchange-correlation functional. Compared 
to PW91 27 and RPBE 28 , we found negligible differences 
only. 

The charge distributions and transfers were analyzed 
by means of the Bader analysis.^ 



III. GROUND STATE PROPERTIES OF 

Aui_ 4 /GRAPHENE 

There are different ways to obtain clusters of a certain 
size on a substrate. Here we consider two possibilities: 
first, mass selected clusters soft-landed onto graphene 
and, second, the clustering of smaller fragments, already 
pre-adsorbed on the graphene sheet. 

When exploring the first scenario, we obtain the 
ground state structures of Aui_4/graphene, see Fig. [I] 
(a)-(d), by full y r elaxing the gas phase cluster ground 
state struct ureeP^H i n various ways on graphene. As 
will be shown in Sec. |IV| the same structures for AU2-4 
are formed when starting from smaller pre-adsorbed frag- 
ments, i.e. by the second scenario. 

The tetramer shows an interesting peculiarity. In 
the gas phase the diamond shaped isomer, Auf*, is the 
ground stateP^^l According to our calculations, it is 
43 meV lower in energy than the Y shaped isomer, 
AuJ\ cf. Refs |30|32j . On graphene the situation is re- 
versed: Au^f /graphene is 137 meV lower in energy than 
Auf /graphene. The overall geometry of both isomers is 
conserved on graphene, cf. Fig. [l](d) and (e). In addition 
to soft-landing, both isomers can be formed on graphene 
by clustering of smaller pre-adsorbed fragments, i. e. A ui 
+ Au 3 -> Auf and 2 • Au 2 -» Au^f, see Sec. |Tv| for 



more details. Note that two of the bonds in A1I4 are 
comparable to the inter-dimer bond, i.e. 2.55 and 2.6 A 
compared to 2.67 and 2.72 A, cf. Figs[l](b) and (e). 



A. Binding sites and energies, charge redistribution 

The adatom, dimer, and Y shaped tetramer prefer the 
binding site on top of a carbon atom, whereas the trimer 
and diamond shaped tetramer bind right above a carbon- 
carbon bond, see Fig. [I] for an illustration. It is inter- 
esting to observe a trend to minimize the number of gold 
atoms binding to graphene for all five clusters, i.e. only 
one Au atom is close to the graphene layer. This al- 
ready indicates that the Au-Au interaction is stronger 
than Au-C interaction, which will be discussed in greater 
detail below. 

The ground state structures of Aui 5 2 /graphene, see 
Fig. [I] (a) and (b), resemble the adsorption behavior 
of Aui 5 2 on MgO and graphite P 8 | 33 | S4 1 The adsorption 
behavior of the trimer and tetramers on graphene, Figs 
[I] (c)-(e), is in contrast to the adsorption of small gold 
clusters on MgO terraces. On MgO the gold trimer and 
tetramer form two bonds between gold and surface oxy- 
gen. Furthermore only the diamond shaped AU4 isomer 
is stable on that substrate.^ 

Care must be taken when comparing our results with 
earlier studies of Aui_ 4 on HOPGpJ^S] First of all the 
authors relied on the use of the local density approxima- 
tion (LDA), which only coincidently describes the bulk 
properties of graphite better than GGA. Secondly, the 
number of k-points was eight at best, making a meaning- 
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ful calculation of, for instance, the DOS impossible. Last 
but not least, spin-polarization was neglected in Refs 
|18|20| . We find spin-polarization especially important 
for open shell systems like Aui^/graphene. They have 
two distinct peaks, for spin up and down states, respec- 
tively, close to the Fermi level, cf. Figs [4] and [5] Albeit 
all these methodological differences one can conclude that 
the ground state structures of Aui_3 adsorbed on HOPG 
and graphene seem to by very similar. Note that Au3 ; 4 
have not been studied on graphene before. The second 
isomer, Au^f Fig. [l](e), has not even been considered on 
graphite before. Note that the tetrahedral structure of 
AU4 is known to be unstable in the gas phase^MI anc [ 
neither did we find this structure on graphene. 




-4.00 



FIG. 1: (Color online) Ground state structures of Aui_4 on 
a 5x5 graphene sheet. Structures (a)-(d) were obtained by 
soft-landing the clusters from the gas phase. The isomer AuJ 
(e) is formed by clustering two pre- adsorbed dimers. It is 
137 meV lower in energy than Auf /graphene (d). Also given 
are relevant bond lengths and Bader charges. Note: in our 
calculations neutral Au atoms have a Bader charge of -11 and 
C atoms of -4, respectively. 

Figure [5] (a) shows the cluster adsorption energies i? a ds 



of the Aui_3 ground state structures, as well as of both 
AU4 isomers. Au^f has the highest E a( ± s of -0.59 eV and 
the adatom having with -0.1 eV the lowest. As mentioned 
above, Au^f is not only the ground state structure on 
graphene, but also adsorbs stronger on it than Auf*, i.e. 
^ads = -0.59 vs. -0.41 eV, respectively. 
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FIG. 2: (Color online) Upper panel (a): Cluster adsorption 
energy of Aui_4 on graphene. Lower panel (b): total charge 
transferred from the graphene sheet into the clusters, calcu- 
lated by the Bader analysis. Full symbols: gas phase struc- 
tures adsorbed on graphene, cf. Figs [I] fa) - (d). Open sym- 
bols: Y shaped Auf /graphene, cf. Fig. |l|(e). 

In all the five cases the whole gold cluster receives, 
according to our Bader analysis, a small charge of up 
to -0.12 e~ from graphene, see Fig. [2] (b). Hence ad- 
sorbing Aui_4 on graphene corresponds to a p-doping of 
the material. A clear odd-even oscillation, depending on 
the number of gold atoms in the cluster, can be seen in 
the charge transfer. It is generally known that the elec- 
tronic characteristics of low-dimensional gold structures 
can exhibit odd/even oscillations depending on the par- 
ity of the number of Au atoms, w hich are determined by 
the opening /closing of the s-shell F ' 36 ' 37 ' 

In addition to the charge transfer into the clus- 
ters, there occurs also a significant charge redistribution 
within the clusters themselves, see Bader charges of the 
individual atoms in Fig. [T] The insets in Figs [4] - [8] also 
illustrate it by showing isosurfaces of the charge density 
redistribution, i.e. Ap = p(Au n /graphene) — p(Au n ) — 
p(graphene), upon the cluster adsorption on graphene. In 
the case of AU2-4 up to 0.19 electrons are redistributed 
from the gold atom binding to graphene into the rest of 
the cluster. 

Concerning the adsorption strength of these small gold 
clusters on graphene we can summarize that two effects 
exist, which can compete with each other. First, the at- 
traction due to the formation of chemical bonds between 
gold and carbon. Second, there also exists a Coulomb in- 
teraction between the charged gold atoms in the cluster 
and the carbon atoms in the proximity of the adsorp- 
tion site. Mostly this Coulomb interaction is repulsive 
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in nature, as it can be seen from the extra charge in the 
clusters that tends to be concentrated as far away from 
the graphene sheet as possible, see Figs [I] (a) - (c). 

Although the extra charge density concentrates on the 
two outer atoms in Auf* /graphene, those two atoms are 
repelled away from the graphene sheet, recognizable from 
the 5% stretched bonds with the bottom gold atom, 
see Fig. [I] (d). In contrast to these four systems, the 
Coulomb interaction between the adsorbed Au^f and the 
underlying graphene sheet is predominately attractive, 
i.e. between the left-most gold atom, cf. Fig. [l](e), car- 
rying most of the additional charge, and the positively 
charged carbon atoms under it, cf. inset in Fig. [7[ 

The comparatively small adsorption energy of the sin- 
gle gold atom is also explained by these competing effects, 
see Fig. [I] (a), which is manifested in the 5-9 % longer 
Au-C distance of Aui /graphene compared to the other 
clusters. 
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FIG. 3: (Color online) Buckling, i.e. vertical distortion, of 
the first three coordination shells around the adsorption sites 
of Aui-4. Full symbols: gas phase structures adsorbed on 
graphene, cf. Figs[l](a) - (d). Open symbols: Auf /graphene, 
cf- Fig. 0(e). 



It has, predominantly qualitatively, been reported 
in the l iterat ure that the adsorption of Aui_4 on 

HOPcpmnnn and of Aui 2 on graphem i39i4oi leads to a 

shift of the carbon atoms in the vicinity of the adsorption 
site towards the gold. Due to the strong sp 2 interatomic 
C-C bonds the vertical distortion is not solely restricted 
to the carbon atoms underneath gold. Instead it spreads 
at least to the third coordination shell around the adsor- 
bate, cf. Fig. [3j Lets consider for the moment only the 
gas phase structures adsorbed on graphene. The full sym- 
bols in Fig. [3] show that the vertical distortion of the car- 
bon atoms binding to the Au atom scales approximately 
linear with the size of the gold cluster. On the other 
hand, the buckling in the second and third coordination 
shell does show a clear odd-even oscillation depending on 
the number of Au atoms, Fig. [3] (full symbols). But the 
actual tetramer ground state geometry Au^f /graphene 
breaks these trends, showing a significantly smaller dis- 



tortion of the C atom positions than Au^ /graphene, see 
Fig. [3] (open symbols). 

Along with the buckling comes a polarization of the 
graphene sheet, i.e. charge deviations from otherwise 
neutral charge state of the carbon atoms, of up to ± 
0.18 e~ . For the systems under consideration the polar- 
ization do es no t show such a clear pattern as observed 
on HOPGP^oj p or instance, the carbon atom binding 
to Aui ? 3 an d Auf gain additional electrons, whereas it 
stays neutral when binding to Au^f or even looses some 
charge as under Au2- Also in the second and third coor- 
dination shells around the binding sites a general pattern 
could not be observed, i.e. charge losses and gains could 
be observed within the same coordination shell. 



B. Electronic structure of Aui^/graphene 

Both the gold adatom and trimer possess an unpaired 
6s electron that gives rise to a total spin moment of 
in the system. 




Au/graphene spin up (down) 

graphene x 0.5 
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FIG. 4: (Color online) Projected spin-polarized total DOS of 
Aui /graphene (black lines; solid spin up; dashed spin down) 
and for comparison pure graphene (red line). Note: the Fermi 
energies of both systems were shifted to eV. Inset: charge 
density redistribution due to adsorption of Am on graphene; 
increase (decrease) of the charge density in red (blue). The 
shown isosurfaces correspond to ± 7-10 -4 e~/A 
on the top show the atomic energy levels of 
atoms. 



3 The arrows 
single gold 



Figures [4] and [5] show the total spin-polarized DOS 
of Aui 5 3 /graphene and, for comparison, also the total 
DOS of pure graphene. For illustration the atomic energy 
levels of an isolated gold atom are shown as arrows in Fig. 
[4] Only in this simple case the origin of the electronic 
states of the adsorbed adatom can easily be identified 
from the gas phase species. 

From the Im and site decomposed DOS (not shown) 
we identified the character of the labeled structures 1-4 in 
the DOS of Aui /graphene and accordingly for the bigger 
clusters as well. 
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In Fig. [4] the structures no. 1 (2) are gold 5<i-states 
with spin polarization up (down). The higher peaks in 
1 and 2 consist partially of d z i -states that overlap with 
p^-states of the underlying carbon atom. Also the gold 
6s-state near the Fermi level, i.e. peak 3, hybridizes with 
the p^-state to form a a-bond. The unoccupied states 
in peak 4 are of similar character as those in no. 3. 
Compared to the atomic states, i.e. arrows 1' and 2' in 
Fig. [4j the d-states are split upon adsorption of the gold 
atom on graphene. Furthermore they are shifted approx- 
imately 0.7 eV downwards in energy. The gap between 
the highest occupied states and the lowest unoccupied 
states AEg, i.e. between the peaks 3 and 4, is 0.15 eV 
narrower than in the gas phase, i.e. 3' and 4'. 

For the bigger clusters, we concentrate our efforts on 
the identification of the most important features of the 
densities of states, i.e. Figs[5]-[8j and their contribution 
to the formation of the clusters and their bonding to 
graphene. 

Already at a first glance at the DOS of AU3 /graphene, 
Fig. [5j one can recognize the higher complexity of the 
structure. The two peaks labeled with 1 comprise pre- 
dominantly of d yz and d x i _ y i -states of the gold atom 
closest to carbon, hybridizing with its p-states. The 
peaks no. 2 are predominantly d x i _ y 2 -states located at 
the two gold atoms at the top. The main features, no. 
3, mostly contain (i-states of all characters from the top 
gold atoms. Most of the d z i -states of all three gold atoms 
are concentrated in this peak, forming the chemical bond 
to carbon by hybridizing with its p^-states. Peak no. 4 
consists of intracluster d xz -d z i states. The peak no. 5 
directly below the Fermi energy contains the 6s-states 
of the top gold atoms that interact only indirectly, via 
(i^-states of the gold atom closest to graphene, with the 
carbon p^-states. Clearly, most of the spin-moment of 
the whole system is concentrated in these states. Finally, 
the peaks 6 and 7 in the unoccupied DOS are overlaps of 
empty 6s-states and carbon p^-states from the top gold 
atoms and the bottom one, respectively. 

From the shown DOS, especially the peaks near the 
Fermi energy, one can conclude that even a small external 
bias of less than ± 0.5 V will significantly increase the 
electrical conductivity of these systems. 

The charge density redistributions upon adsorption Ap 
are shown as insets in Figs[4]and[5j The earlier mentioned 
polarization of the carbon atoms in the first three coor- 
dination shells around the adsorption sites can be seen. 
The extra charge in Aui also leads to a charge density 
redistribution within the gold atom, whereas the gold 
atoms in the trimer seem to equally share the additional 
charge from the substrate. Therefore it is important to 
take the Bader analysis, Fig. [I] into consideration as well 
that showed a charge depletion of the bottom gold atom. 
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FIG. 5: (Color online) Projected spin-polarized total DOS of 
Aus/graphene (black lines; solid spin up; dashed spin down) 
and for comparison pure graphene (red line). Note: the Fermi 
energies of both systems were shifted to eV. Inset: charge 
density redistribution due to adsorption of AU3 on graphene; 
increase (decrease) of the charge density in red (blue). The 
shown isosurfaces correspond to zb 7T0 -4 e~/A 3 . 



C. Electronic structure of Au2,4/graphene 

Obviously, the total magnetic moment of the closed 
shell structures Au2,4 must be zero in the gas phase. This 
is unchanged by adsorbing them on graphene. 




-8 -6 -4 -2 
Energy (eV) 

FIG. 6: (Color online) Projected total DOS of Au2/graphene 
(black line) and for comparison pure graphene (red line). 
Note: the Fermi energies of both systems were shifted to 
eV. Inset: charge density redistribution due to adsorption of 
Au2 on graphene; increase (decrease) of the charge density 
in red (blue). The shown isosurfaces correspond to db 7T0 -4 
e-/A 3 . 

Upon adsorption the symmetry of the electronic struc- 
ture of the dimer is broken, see Fig. [6j From the Im 
and site decomposed DOS (not shown), we see that peak 
no. 1 consists of d z i -states. The broader feature no. 2 
and peak no. 3 are d xy , d yzi and d x i_ y i -states, predomi- 
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nantly from the gold atom closest to carbon, hybridizing 
with carbon p XyV -states. Peaks 4 and 5 do not contribute 
to the chemical bonding to graphene, as they are s-d- 
states localized at the cluster. Whereas no. 6, mostly s- 
d z 2 -states at the top gold atom, indirectly interacts with 
carbon p z via d z i at the bottom gold atom. Finally, 
the two peaks at no. 7 are predominantly unoccupied 
6s-states of the gold atom close to carbon, overlapping 
empty p^-states of carbon. Note also that the states in 
peaks 4-7 substantially narrow the cone in graphene's 
DOS around the Fermi energy, leaving it only undis- 
turbed in the interval Ep ±0.5 eV, Fig. [6j 



T 




Energy (eV) 

FIG. 7: (Color online) Projected total DOS of the Y shaped 
Au4/graphene (black line), see Fig. [I] (e), and for compar- 
ison pure graphene (red line). Note: the Fermi energies of 
both systems were shifted to eV. Inset: charge density re- 
distribution due to adsorption of Au4 on graphene; increase 
(decrease) of the charge density in red (blue). The shown 
isosurfaces correspond to =b 7-10 -4 e~/A 3 . 

To summarize the predominant nature of the labeled 
peaks in the DOS of Au^f /graphene, Fig. Rwe label 
the gold atoms in the following way, see Figr"[l](e): Au^ 
(left), Au m (middle), Au £ (top), and Au 6 (bottom). The 
carbon atom under Au 6 , forming the chemical bond with 
the cluster, is named C 6 . 

Peaks no. 1 and 2 consist of s-d-states, d z i -states re- 
spectively, localized on the cluster only. Peak no. 3 con- 
sists of d-states that contribute to the chemical bond by 
overlapping with p y -states of C 6 , i.e. -states on Au m 
and d^-states on Au^' 6 . The peak no. 4 and its shoul- 
der consist mainly of d xy and d x i_ y i on Au 6 overlapping 
with carbon p-states as well. Peak no. 5 are d x 2_ y 2- 
states solely localized within the cluster, while no. 6 are 
s-d z 2 -states of Au 6 hybridizing strongly with p z of C 6 . 
The highest peak no. 7 consists of d xy and d x i _ y 2 -states 
of Au^, Au m , and Au £ , i.e. a major part of the strong 
gold-gold bonds that hold the cluster together. Peaks no. 
8 and 9 are d-states almost exclusively localized on Au m 
and Au^, respectively, as well as peaks no. 10 (d x 2_ y 2- 
states) and 11 (s-d^-states) which are localized on Au^. 
As we saw for the other clusters the unoccupied states at 



peak no. 12 are s-states, in the present case of Au £ and 
Au 6 , overlapping with unoccupied p^-states of C 6 . 
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FIG. 8: (Color online) Projected total DOS of the diamond 
shaped Au4/graphene (black line), see Fig. [I] (e), and for 
comparison pure graphene (red line). Note: the Fermi en- 
ergies of both systems were shifted to eV. Inset: charge 
density redistribution due to adsorption of Au4 on graphene; 
increase (decrease) of the charge density in red (blue). The 
shown isosurfaces correspond to =b 7T0 -4 e~/A 3 . 

Despite its more symmetric geometry the DOS of 
Au^ /graphene, Fig. |8j has many features which can be 
recognized from the other isomer, cf. Fig. [7| For conve- 
nience we label the gold atoms in Auf /graphene again: 
Au 6 (bottom), Au m (middle), and Au £ (top), see Fig. 
[I] (d). The two carbon atoms, forming bonds with the 
cluster, have a very similar electronic structure and we 
therefore simply name them both C 6 . 

The relatively broad peak no. 1 consists mainly of d xz - 
states of Au 6 and Au £ . The shoulder below peak no. 2, 
are d^-d^-cluster states of Au 6 and Au £ , with a smaller 
amount of intermixed s-states that hybridize with the p- 
states of C 6 . The actual peak no. 2 and the double-peak 
no. 3 are predominantly d yz and d x 2_ y 2 -states of Au 6 
and Au^, respectively. The lower peak no. 4 contains 
s-d xy -sta,tes of Au 6 that also overlap with C 6 p^-states. 
The high peaks no. 5 to 7 are the strong intracluster 
bonds, formed out of d-states. Peak no. 8 on the other 
hand consists mainly of d^^-states of Au m , as well as 
d z 2 -states of Au^ and Au 6 , contributing to the chemical 
bond to carbon by hybridizing with p^-states of C 6 . The 
highest occupied states, peak no. 9, are again intracluster 
bonds, formed out of s-states of the Au m and cL^-states 
of Au^ and Au 6 . As for the other clusters the lowest 
unoccupied states, peaks no. 10, are unoccupied s-states 
on all the four gold atoms, that overlap with unoccupied 
C 6 p 2 -states. 

The gap between the highest occupied states and the 
lowest unoccupied cluster states AEq in Fig. [7] and [8] 
is very similar, but smaller than that in the dimer case, 
Fig. [6] Hence, in contrast to Aui^ /graphene, an external 
bias of more than ±0.5 V is needed to see a significant 
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increase in the electrical conductivity. 

IV. MOBILITY AND CLUSTERING 

In this section we discuss the mobility and initial clus- 
tering processes of gold on graphene. First experiments 
on the diffusion and desorption of gold atoms on graphite 
date almost three decades back. 15 With the rise of STM's 
quantitative estimates of the diffusion and desorption 
barriers became feasible. These experiments revealed 
that the binding energies of Aui^ on graphite exceed 
their diffusion barriers!^ Computations of gold adatoms 
and dimers on HOPG supported these experiments. 19 



stable binding site. These transition states define the dif- 
fusion barriers of the gold clusters along the C-C bonds. 
For convenience the diffusion barriers are also explicitly 
given in Fig. [9] They range from only 4 meV for Au2 
and Auf^, making them highly mobile, to 36 meV for the 
least mobile trimer. 

The adsorption energies of all five gold clusters on 
graphene, cf. Fig. [2] (a), exceed their diffusion barriers 
by at least one order of magnitude. Hence, while these 
small gold clusters will rather strongly adsorb onto the 
graphene sheet, they will easily diffuse even at very low 
temperatures. 



B. Clustering of pre-adsorbed fragments 



A. Diffusion along the C-C bonds 

We determined the diffusion barriers for Aui_4 on 
graphene from the total energies on the top (t) and bridge 
(b) positions, as well as from one position in between 
those two (tb). These energies are sufficient, since none 
of the five clusters was found to bind to graphene at the 
hollow sites, i.e. in the center of a carbon hexagon. 
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FIG. 9: (Color online) Difference in the total energies AEq 
at the top (t), bridge (b), and a positions in between (tb) of 
Aui-3 and Au 4 ' on graphene. Eo of the respective cluster 
ground state was subtracted to show the diffusion barriers 
along the C-C bonds (also given in parenthesis). Note the 
smaller scales for Au3,4- 

Figure [9] shows the differences in the total energies, 
AEo, of the five clusters with respect to the total energy 
at their binding site, i.e. top (Aui 5 2 and AuJ) and bridge 
for the remaining two. 

Note that only Au2 binds on both the top and bridge 
position, while the other four clusters have a saddle point 
at the bridge (Aui and AuJ) or top (AU3 and Auf) 
position, respectively, i.e. if allowed they relax to their 



We modeled the initial clustering processes of gold on 
graphene in two steps. First, we simulated the formation 
of a dimer out of two pre-adsorbed gold atoms, cf. Figs 



10] and [TT] Second, two possible pathways for forming a 
tetramer, cf. Fig. [l2j were explored. At the same time, 
the formation of a trimer out of a dimer and adatom was 
studied, too. 

In order to simulate the clustering of two pre-adsorbed 
gold atom, they were initially placed as far away from 
each other as possible on a 5x5 graphene sheet, i.e. 7.1 A. 
This configuration is named position 1 in Figs [To] and 
[TT] Since the total energy of the system decreases only 
marginally between position 1 and 2, we omitted inter- 
mediate steps, cf. Fig. [To] 

In the subsequent steps, i.e. position 2 to 6, the sec- 
ond gold atom was brought closer to the fixed one along 
two C-C bonds, see inset in Fig. [Tol For each of the 



positions 1-6 only the x-y-coordinates of the adatoms 
and one carbon atom at the rim were kept fixed, while z- 
coordinates of the adatoms and the rest of the graphene 
sheet could fully relax. In a final step we allowed the 
system to freely relax from position 6 and without an 
additional activation barrier it converged to the dimer 
ground state (position 7). This final structure is identi- 
cal to the one we found in Sec. |III| by soft-landing the 
dimer on graphene, see Fig. [l](b). 

Both, in the gas phase and on graphene, the formation 
of a gold dimer is exothermic with a gain in total energy 
per gold atom of AEq = -1.14 and -1.29 eV, respectively, 
see Fig. 
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The higher energy gain on graphene is a 
result of the higher adsorption energy of Au2 compared 
to two adatoms, see E a a s in Fig. [2] (a). AEo largely ex- 
ceeds both the adsorption energy as well as the diffusion 
barriers of Aui /graphene. Thus the driving forces be- 
hind the clustering of two pre-adsorbed gold adatoms is 
the opportunity to lower the total energy of the system 
by forming strong gold-gold bonds, i.e. the strongly hy- 
bridized d-states in Fig. [6] At the positions 3 to 6 forces 
act on the two adatoms that try to pull them towards 
each other. This attractive, graphene mediated interac- 
tion has a maximal range of approximately 4.3 A, i.e. the 
Au-Au distance in position 3. At this distance the two 
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gold atoms are still well separated in the gas phase. 

Note that the last two positions represent a dimer that 
is parallel (position 6) and perpendicular (position 7) to 
the graphene sheet. Obviously, in the gas phase the total 
energies have to be identical for position 6 and 7. 




1 2 3 4 5 6 7 
Positions of the Au atoms 

FIG. 10: (Color online) Changes in the total energy AEq 
when forming Au2 /graphene from two pre-adsorbed adatoms 
(solid symbols) and for the dimer formation in the gas phase 
(open symbols) . The x-y coordinates of the Au adatoms were 
kept fixed for positions 1-6. The system relaxes without an 
additional activation barrier from position 6 to 7, the dimer 
ground state structure (Fig. [5](b)). The Au-Au distances for 
the different positions are given on the top axis. 

Figure [TT] shows the changes in the spin-projected to- 
tal DOS during the intermediate steps of forming a gold 
dimer on graphene, cf. position 1 to 7 in Fig. [10] 
Note that even for the largest possible separation of the 
two adatoms in the 5x5 super-cell, i.e. 7.1 A (position 
1), there is still a electrostatic interaction between the 
charged gold atoms, cf. Fig. [I] (a). This can be seen 
from their 6s-states near the Fermi level that should ex- 
actly coincide without any mutual interaction. At a Au- 
Au distance of 4.9 A (position 2) the 6s peak splitting 
becomes more visible. But still the 6s-states below the 
Fermi level belong solely to the spin-up channel. When 
the graphene mediated attraction between the adatoms 
becomes obvious at a distance of 4.3 A, see AE at posi- 
tion 3 in Fig. [lOj one of the 6s spins flips, leaving the total 
system with no total magnetic moments, as we saw above 
for the closed shell system Au2/graphene. Upon further 
approach of the two gold atoms the two atomic 6s-states 
form a molecular orbital that moves to lower energies. 
Finally, it hybridizes with the rf^ -states, cf. peaks 4 and 
6 in Fig. [6j in the dimer ground state structure (position 
7). Meanwhile the unoccupied 6s-states directly above 
the Fermi level in pos. 3-6 move up in energy to finally 
form peak no. 7 in Fig. [6] 

Our next aim was to study the formation of Au3 5 4 on 
graphene. We accomplished this by allowing two pre- 
adsorbed fragments to fully relax from an initial state 
corresponding to position 6 in Fig. 10 There are two 
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FIG. 11: Changes in total spin-polarized DOS (solid lines spin 
up; dashed lines spin down) when two pre-adsorbed adatoms 
come stepwise closer to each other to form Au2 /graphene, cf. 
Fig. |10| The diagonal line marks the Fermi energies at the 
different positions. 



different reaction pathways for forming a tetramer, re- 
sulting in the two isomers Auf* and Au^f, see Figsjl] (d) 
and (e). In the first reaction pathway two dimers cluster 
from their perpendicular ground state into the Y shaped 
tetramer, 2 • Au2 —> Au|\ The energy gain per gold 
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atom during this step is AEq — -0.27 eV, see Fig 
(a). Nonetheless, per dimer this energy gain is still big- 
ger than the dimer adsorption energy of -0.47 eV, see Fig. 
[2] (a), overcompensating the breaking of one Au-C bond 
to form Au^f /graphene. 
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FIG. 12: (Color online) Changes in the total energy AEq 
during the formation of Auf /graphene out of (a) two pre- 
adsorbed dimers, and (b) Auf out of four adatoms (solid 
squares). For comparison AE§ for the formation of Auf* in 
the gas phase is shown as well (open circles). The dotted 
horizontal line marks the highest gain in energy. 

The second pathway includes two intermediate steps to 
form Auf* out of four single gold atoms. To begin with, 
a trimer is formed, Aui + Au2 —> AU3. Then, by adding 
yet an adatom, one forms the tetramer: Aui + AU3 — > 
Auf. The energy gain AEq during the whole reaction is 
shown in Fig. 12 (b). This process is exothermic and the 
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fragments only have to overcome their diffusion barriers, 
Fig. [9] As for the dimer formation, the energy released 
by forming the strong Au-Au bonds is the driving force 
behind the clustering, which is for all steps bigger than 
the cost of breaking a Au-C bond or overcoming diffusion 
barriers. 

Although we did not explicitly calculate the stepwise 
movement of the different fragments, as we did in the 
case of two adatoms, see Fig. [TUJ we expect a similar 
behavior due to the high mobility of the pre-adsorbed 
clusters, see Sec. |IV| 

A comparison of Aifo in Figs [121(a) and (b) shows, as 
mentioned in Sec. Ill, that Auf/graphene is 137 meV 
lower in energy than Au^/graphene. The open symbols 
in Fig. [12] show the energies of the reaction steps in 
the gas phase. Since Auf is the gas phase ground state 
structure, we assumed that both reaction pathways give 
this structure. The end points in Fig. 12 (b) seem to 
coincide, i.e. the energy gain by clustering in the gas 
phase and on graphene appears to be the same. That is 
actually the case since the adsorption energy of four Aui 
on graphene is -0.4 eV, which is almost identical to E a( ^ s 
of Auf with -0.41 eV. 



V. SUMMARY AND CONCLUSIONS 

By means of density functional theory, the adsorption 
of Aui_4 clusters on graphene, their mobility and clus- 
tering from pre-adsorbed fragments have been studied. 
This allowed us to explain the initial steps in the exper- 
imentally observed formation of bigger gold aggregates 
on graphene. While Aui_3 adsorb similarly on graphene 
as on graphite, the tetramer turns out to be more pecu- 
liar. In the gas phase it has two isomers, the diamond 
shaped Auf^, being the ground state, and the Y shaped 
Au|\ On graphene Au^f becomes the ground state in- 
stead. Although Au^f binds stronger to the substrate 
and also receives more charge from it than Auf does, it 
disturbs the graphene sheet substantially less. 

The cluster adsorption energies on graphene of all stud- 
ied clusters range from -0.1 to -0.59 eV and substantially 
exceed their diffusion barriers, which are 4 to 36 meV, 



only. All clusters diffuse along the C-C bonds, since none 
of them binds to the hollow sites in the centers of the car- 
bon hexagons. 

Our detailed analysis of the densities of states shows 
which states contribute to the strong intracluster bonds, 
predominantly the 5d- states, and which to the chemical 
bond between gold and carbon, i.e. mostly s- d z 2 -cluster 
states hybridizing with p^-states of carbon. 

The formation of the strong Au-Au bonds is the driv- 
ing force behind the tendency to form bigger clusters 
from smaller pre-adsorbed fragments. We studied these 
processes in detail for the stepwise formation of a dimer 
out of two pre-adsorbed adatoms, including the evolu- 
tion of the electronic structure during the clustering. The 
exothermic reaction has no additional activation barriers 
apart from the small diffusion barriers of the fragments 
moving along the C-C bonds. We find the graphene me- 
diated, attractive interaction of two gold adatoms to have 
a range of approximately 4.3 A, which is clearly longer 
than in the gas phase. 

Finally, we have also studied the formation of the 
trimer and tetramers by the reactions Aui + Au2 —> AU3, 
Aui + AU3 — » Auf*, and 2 • Au2 — » Au|\ Both reactions 
are also activation barrier free. As for the dimer, the 
energy gained by clustering, i.e. forming Au-Au bonds, 
exceeds substantially both the Au-C bond energies as 
well as the diffusion barriers. Note that both tetramer 
isomers, i.e. Au% ,D are formed on graphene, depending 
on the actual reaction pathway. 
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